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Abstract. We investigate the analytically solvable pion-laser model, and its gen- 
eralization to arbitrary multiplicity distributions. Although this kind of extension 
of the model is possible, the pion laser model in its original form is unique: it is 
the only model in its class that posesses an analytic solution. 

This work is dedicated to the memory of the late J. Zimdnyi and corresponds to a follow-up 
work on one of his favourite topics: the analytic solution of the pion laser model. 



1 Introduction 

Correlation between different identical particles from a thcrmalized source is due to the bosonic 
or fermionic nature of them. The effect was discovered by R. Hanbury Brown and Twiss in 1954, 
they pointed out that the observed correlation of photons carry information on the angular 
diameter of the emitting distant star fll. In experimental high-energy physics, the phenomenon 
is known as GGLP effect [2 , one commonly uses the term HBT-effect for the fact, that boson 
correlations appear and the shape of the correlation function is related to the geometry of the 
source. Since its discovery the HBT effect was proven to be an important tool in high energy 
physics for investigating the space-time extent of the reaction process. 

In heavy ion physics, experiments usually measure charged pion correlations. If one ne- 
glects the final state interaction (FSI) and multi-particle correlation effects, the description of 
the correlation is simple: the correlation function is essentially the Fourier transform of the 
emission source distribution, referred to as source function. But in realistic situations neither 
FSI nor multi-particle correlations are negligible. Two-particle FSI can be taken into account 
either by generalized Coulomb corrections or by using the so-called imaging method [3] . Multi- 
particle correlation means that the final state outgoing wave-function needs to be symmetrized 
in all variables. This obviously requires n! operation, where n is the number of particles. If 
the phase-space density of produced pions increases, then the effects arising from multi-particle 
symmetrization become significant, so at first sight this is a non-polynomially (NP) hard prob- 
lem. Note that the computing time of NP hard problems of the order " n" increases faster than 
any power of n. Hence these problems for large n are essentially unsolvable numerically. 

But in fact, there are such models in which one can overcome this difficulty and calculate 
the correlation functions using fully symmetrized wave- functions. The first success in this di- 
rection was achieved by S. Pratt [6j, who showed that in a special case of source functions and 
factorization, the NP-hard problem reduces to a set of recurrence relations, called ring algebra, 
detailed first in ref. 7]. The next step were refs. |4|5j . where the model was generalized to wave- 
packet states instead of plane-wave final states, and the recurrence relations were solved. The 
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model had only one flaw: the multiplicity distribution of it is very special, and in some sense, 
unphysical. For example, coherent behaviour corresponds to Poisson multiplicity distribution 
in case of optical the optical lasers, while fully chaotic, rare gas limit corresponds to the Poisson 
multiplicity distribution in case of the pion laser model. In the present work we investigated if 
the original model of refs. [4151617] can be modified in such a way, that the Poisson multiplicity 
distribution will characterize the pion condensate, and not the chaotic state. 

In sections[2l[3l|4]and[5]we recapitulate and slightly re-formulate earlier results of refs. 4 "5"6'7j 
to prepare their generalization, utilizing the idea of Poisson transformation. Section [S] describes 
the outcome of this investigation and Section [7] summarizes the results. 



2 Definition of the model 

The investigated pion- laser model describes a multi-particle system containing arbitrary number 
of bosonic wave-packet states. The notations ^, x ^-nd stand for the center in coordinate 
space, the center in momentum space and the width in momentum space of a given wave 
packet, respectively, and a :— (^, x, a) is an abbreviated notation for these prarmeters (for 
example, integrating over a means integration over all variables). The quantities ^ and x £^re 
three- vectors. Such a one-particle state |a) characterized by the parameters a = (^,X,cr) is 
created from the |0) vacuum state by d^, the following wave-packet creation operator: 

exp\-^{p-xf -ii.{p~x)\a\p), (1) 



where a\p) is the usual pion creation operator. The normalization of the states \a) is {a\a) — 1, 
while states with different parameters have overlap, which we denote by {ai\aj) :— ^ij. With 
this, the properly normalized n-particle states are 

Crt \ -1/2 n 

cr i=l I i=l 

where the summation runs over all the possible a permutations of the n indices, and cr^ means 
the image of the i index by the permutation a. We note the expression of jij (if cr^ = aj = cr): 

7y = exp <^ — + 2 - Xj) (6 + Q ) ■ (3) 

The investigated model is defined through its density matrix p as 

oo 

/5 = ^ VnPn , Trp = 1 (4) 

n=0 

with pn being the multiplicity distribution, and /5„ the density matrices of events with fixed 
number of particles. The normalization is as usual: Tr/3„ = 1 for any n, so X^^o^*" ^ 
Pn matrices are built up of the multi-particle states defined above: 

j '^'^i • • •da„p(Q:i . . .an) \ai . . .a„)(ai . . .a„|, (5) 

where p {ai . . . an) are density functions. There is a special type of them, in what case the 
model is analytically solvable [4J, namely, if the normalization factor in Eq. ^ is cancelled: 

1 ( " \ " 
Pn{ai...an) = I HIl^J-. I n/'i(")- ('^) 



cr j = l 
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pi(ai) = (5(cri - cr) X — exp<{ -7^1 X exp (7) 

with T being the temperature and i? the radius. The normalization factor in Eq. ([6]) is 

M{n) = ^'W\ ^a^Pl (a^)7^^, • (8) 

a i=l 

The definition of the model is completed with the specification of the multiplicity distribution 
Pn, which step is postponed to Section [6l 



3 Observables 

We want to calculate inclusive n-particle distributions as well as exclusive ones (considering 
events only with fixed m number of particles). The definitions of the exclusive and inclusive 
quantities are, respectively 

TVf") (fci, . . . , fc„) = Tr {p^aHki)a{h) . . . a^k^Mk^)} , (9) 

OO 

7V,(fci,...,fc„) =Tr{pat(fci)a(fci)...at(fc„)a(fc„)} = ^ p„7vf™^ (fci, . . . , fc„) . (10) 

■m—i 

The normalizations are 

d^fc ivf™) (fc, fci, . . . , k,^,) + Ntl (fci, . . . , h-i) , (11) 

d3fci...d3fc, 7V(™)(fci,...,fc,) =r^^.. (12) 

(to — i)! 

d^ki...A^ki N,{ki,...,h) = (n(n~l)...(n-i + l)), (13) 

where ((n) . . .{n — i + 1)) is the i-th factorial moment of the multiplicity distribution. Because 
of the special choice of the density functions in eq. ^ in the present model it is possible to write 
down a relatively simple expression for the exclusive distributions with the auxiliary function 



p{Pi q) = J da {a\p)pi{a){q\a). (14) 
Performing some combinatorics one finds that 

From this one obtains the lower order exclusive functions as 

7Vi")(fci,...,fc„)^ ) X /d3fc„+i...d3fc„vf[p(fc„A:,J. (16) 

[n — my. Al [n] J ^l^±± 

For sake of clarity we write down the form of the two lowest order functions: 

/n 

N^-\k,M)^%r^ /d%...d3fc„vfrp(fc„fc..). (18) 
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The calculation of these n!-termed integrals is usually a NP-hard problem. However, in the 
present model a nice analytic solution is possible by means of a so-called ,,ring algebera" [6]. 
The generalization to the case of wave-packets is found in Ref. [4], in the next section we 
summarize the results. 



4 Analytic solutions 

The integrals in eqs. (fT7|) and (fT5|) can be evaluated using Wick's theorem: the p functions 
in a term of te sum over the permutations can be ordered to rings (i.e. orbits of the given 
permutation) . The evaluation of the integrals then reduces to the evaluation of this so-called 
ring algebra. We just quote the results: one introduces the set of functions Gi {p, q) as 

Gi{p,q) =p{p,q) , Gn{p,q) = J d^Pip{p,Pi)Gn-i{pi,q)- (19) 

These functions appear as ,,ring integrals". Another important quantities are the normalization 
constants Af{n). Following Ref. [J] (but slightly deviating from its conventions), we use an 
alternative notation as ojn = [n). Putting all this together, the one- and two-particle 
exclusive distribution functions are given by 

7Vf")(A;i) = — VG,(fci,fci)o;„_„ (20) 

^ n i—1 

N^''\h,k2) = — ^^{Gi(A:i,A:i)G,_,(fc2,A:2) + Gz(fci,fc2)G,_z(A:2,fci)}c^„-,. (21) 



" i=2 1=1 



The general iVj'"-* could be expressed in a much similar way. The functions G„ {p, q) are found to 
have an analytic form |4j . The key point is the stability property of the Gaussian distribution: 
the wave-packets as well as the density functions in eqs. ([5]) and ([6]) are Gaussians, so the 
G„ {p, q) functions will be also Gaussian, so we parametrize them as 

Giip, q) = h exp {-ai {p^ + q^) + gipq] , (22) 

and then from the recursive definition in eq. p9p one can derive the recurrence relations, which 
determine the quantities /i^, and gi. 

hn+l^hihn( ) , /ll = — TTTi (23) 



4 (ai -I- a„) ' 2(7, 



9n+l = 777 ^ '9l=Reff- (25) 

2(ai+a„) 

The constants ctt and Ref / in the initial values can be considered as effective values of the 
source parameters: 

crl=a^ + 2mT , Rj^^ = + (26) 



T 



The solution of the recurrence relations is given by the following equations: 



9iYn+i , ,„ r/2^\" 1 



3/2 



9n = 1^, (27) 
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where the auxiUary quantities are 



K = ^— , r± = ^±Jf^) -1. (28) 



91 V V 51 

These quantities give the auxiliary d functions. We do not go into the details of the derivation, 
just recommend Ref. [4] for a thorough discussion. We should note that the definitions and 
equations used in the present paper differ slightly from those of Refs. |4|5j . as the multiplicity 
distribution is introduced here in the last possible step. 

For the exclusive spectra one needs the ujn quantities. They obey the following relation, 
obtained from eq. (|20p using the normalization condition of A'^i (k): 

1 " f 

(jJn = -/ l^n-iiCi , C'n = / dfcG„ (/c, fc) . (29) 

n ^ — ' / 

i=\ ■' 

The solution to this can be obtained from the expansion of the following power series: 

OO / CXD \ 

~g{z) = exp C„z" . (30) 

ji=0 \ji=l / 

We do not go in details now, just mention that in two limiting cases the w„ quantities have a 
simple form: if R^^j:(j^ oo (that is, in the rare gas limit) 

C„ = Sn,l =^ UJn = (31) 

n! 

while in the opposite, dense gas limiting case (that is, R^^^a^ 0) 

C„-i ^ c^„ = l. (32) 
n 



5 Inclusive spectra 

In a real experimental situation one usually measures the inclusive spectra, exclusive measure- 
ments can be done much harder. However, for the inclusive spectra ne needs the exclusive ones, 
and so one needs to specify the multiplicity distribution, and at this point difficulties occur. 

We could specify any kind of multiplicity distribution, but eq. (|30)) suggests a way that 
enables us to continue the calculation analytically. Namely, the role of the C„-s in eq. ([29|l is 
very similar than that of the so called combinants of a multiplicity distribution: in this sense 
LOn corresponds to the multiplicity distribution itself. So one reasonable possibility is to choose 
Pn as 



• oo 

E 



(33) 



where no is some ,,mean" value. We note two important relations concerning infinite sums: 

(34) 




OO \ OO / OO \ 

(=1 / fc'=i \i'=i / 



k'=2 




-I'h'-q'Cq' . (35) 
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Using these, the inclusive spectra are found to have a simple form [4] . Introducing the function 

oo 

G{p,q) = Y,Gn{p,q)n^, (36) 



one has 



N,{ki) = G{h,ki), (37) 
N2 (fci, fe) = G iki,ki) G (fc2, fe) + G (fci, fca) G (fcz, ^i) . (38) 

Higher order functions can be obtained in a much similar way. It is important to note two 
specialities of these expressions. First, since G {p, q) has an analytic form, all the measurable 
inclusive distribution functions can be calculated in an easy way. Second, evaluating G (p, q) 
at the intercept point p = q, we get G {p,p) = 1, so, for example, if we define the two-particle 
correlation function as 

r-f \ N2{p,q) ,„„^ 

N,ip)N,iqy 

we see, that it has a physical intercept behavior, C {p,p) — 2. Similarly, for the higher-order 
correlation functions, C„(p, ...,p) = n!. 

This multiplicity distribution was considered in Refs. |4l5j . Our aim was to generalize the 
multiplicity distribution and find analytic expressions for the inclusive spectra. A motivation 
arises from the limiting cases of the present multiplicity distribution: we see from eqs. (|3ip and 
P2p that in the rare gas limiting case p„ approaches a Poissonian distribution with uq mean, 
and in the opposite, dense gas case we have 

We would be interested in a distribution which has the opposite behavior: in the rare gas limit 
it should be a thermal, negative binomial distribution, and in the dense limit a condensed, 
laser-like Poissonian distribution. We will see that this goal can be achieved without difficulties 
by redefining p„, but the analytic simplicity will be lost. 



6 General p„-s, uniqueness 

We can generalize the definition of p„ in eq. (|33p to a case of an integral transformation as 

-1 

^nj/"if(y), (41) 



/ d J V c.„,2/" 

•^0 Vm=0 / 



This is a pretty general multiplicity distribution, since the function iJ (y) is arbitrary. If 
H (y) — 5{y — no), then pn is the original one as in eq. (|33p . In the rare gas limiting case 
the transformation is just 

Pn= / dyexp(-y)^i/(y), (42) 
Jo 

which is the so-called Poissonian transformation of the H (y) function. As an example, let us 
consider the function H (y) = exp (~'{^)' yields a geometrical distribution in the rare 
gas limiting case: 
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Another motivation for the choice of the multiplicity distribution as in eq. (|4ip lies in the 
calculation method of the inclusive spectra. Before going into this, let us mention, that if 
we allow H (y) to depend not only on y, but also on the model parameters {Pefj and cr|,), 
then practically all type of multiplicity distributions can be expressed as an integral transform 
according to eq. (I4ip . So for a given multiplicity distribution one can find the appropriate H (y) 
function. 

We can use eqs. (p4|) and (f35| and their analogues with more variables to calculate the 
inclusive spectra. For example, from eqs. pI7|) and ([M)) A^i (fci) is 

poo TT / \ ^ 

iVi(fci)=/ dy ^ ^;/'^G,(fci,fciV,_, =^G„(fci,fci)/„, (44) 

Jo Z^,„=o^™y ,=1 „=i 

where /„ ~ H (y) y". In the same way, we have for the two-particle inclusive function 

OO OO 

^2 (fci, fe) =J2J2 (^1' (^2, fe) + {h,k2) Gj (fc2, fci)} h+j, (45) 

and so on for Nn {ki, . . . , fc„). 

We can now draw two conclusions 

1. If and only if /„+„ = /„/„, then we can introduce the quantity G (p, q) = J^'iZo Z'*^' (P' 
and in this case the inclusive quantities have a simple analytic form: it is enough to calculate 
G (j), q) once, and all the inclusive functions are simple combinations of it. 

2. The intercept behavior of the correlation functions is physical if and only if fn+m = fnfm- 
That is, for example, the requirement for C2 {p,p) to be 2, is met only if the /„ integrals 
possess the mentioned relation. 

But, fn+m = fnfm holds Only for the case, when H (y) = S {y — no) with some no value. So, 
the conclusion of the present investigation is a negative statement: the multiplicity distribution 
can be generalized, but the analytic simplicity of the inclusive distribution functions as well as 
the physical behavior of the correlation functions do not remain for a more general pn than 
investigated already in Ref. 



7 Summary 

The aim of this analysis was to generalize the pion-laser model. The model was first described 
by S. Pratt [5], and then solved entirely by one of us (T. Cs.) together with J. Zimanyi [4I5| . We 
have found, that although it is possible to modify the definition of the multiplicity distribution 
so that the thermal or chaotic limit corresponds to a negative binomial multiplicity distribution, 
but if we require the analytic solvability of the model, it uniquely prescribes that the multiplicity 
distribiton in the rare gas limit has to be a Poisson. This is the only case, when the intercept 
of the correlation function in the rare gas limit has the expected C2{k,k) — 2 value. In this 
sense, we have proven that the original form of the pion-laser model is unique. A possible other 
path for analytically solvable multi-particle systems could be developed if one omits some of 
the assumptions of the present model: for example, if one factorizes the density function in 
another way than in eq. ([6]). But these generalizations — or more like another models — are 
extremely difficult to handle. 

After the publication of the solution of the pion laser in refs. j4l5l6l7| , a similar generalization 
has been proposed by Heinz, Scotto and Zhang in ref. ^8j. Working with other conventions, they 
also obtained the general form of the inclusive spectra as in eq. (j45|) . however, the emphasis 
of their work was not on the uniqueness of the analytic solution, but more on its numerical 
generalization. With the present work, we in a sense complete the investigations of refs. [4 5j. 
We find that the pion-laser model is unique, in the sense that it is the only the analytically 
solvable model of multi-boson symmetrization of wave-packets in the class of the considered 
models. 
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